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information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite 
1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188,)  Washington,  DC  20503. 

1 .  AGENCY  USE  ONLY  (  Leave  Blank)  2.  REPORT  DATE  30  April  2008 

3.  REPORT  TYPE  AND  DATES  COVERED 

Reprint:  01  July  2007  -  31  June  2008 

4.  TITLE  AND  SUBTITLE 

Space-Time  Power  Schedule  for  Distributed  MIMO  Links  without 
Instantaneous  Channel  State  Information  at  the  Transmitting  Nodes 

5.  FUNDING  NUMBERS 

W911NF0410224 

6.  AUTHOR!  S) 

Yue  Rong,  Yingbo  Hua,  Ananthram  Swami,  and  A.  Lee  Swindlehurst 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  California 

900  University  Ave 

Riverside,  CA  92521 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.  S.  Army  Research  Office 

P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 

10.  SPONSORING  /  MONITORING 

AGENCY  REPORT  NUMBER 

46637. 167-CI-MUR 

11.  SUPPLEMENTARY  NOTES 

The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  be  construed  as  an  official 
Department  of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 

12  a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  Federal  Purpose  Rights. 

12  b.  DISTRIBUTION  CODE 

N/A 

13.  ABSTRACT  (Maximum  200  words) 

A  space-time  optimal  power  schedule  for  multiple  distributed  multiple-input  multiple-output  (MIMO)  links  without 
the  knowledge  of  the  instantaneous  channel  state  information  (CSI)  at  the  transmitting  nodes  is  proposed.  A  readily 
computable  expression  for  the  ergodic  sum  capacity  of  the  MIMO  links  is  derived.  Based  on  this  expression,  which  is  a 
non-convex  function  of  power  allocation  vectors,  a  projected  gradient  algorithm  is  developed  to  optimize  the  power 
allocation.  For  a  symmetric  set  of  MIMO  links  with  independent  identically  distributed  channels,  it  is  observed  that  the 
space-time  optimal  power  schedule  reduces  to  a  uniform  isotropic  power  schedule  when  nominal  interference  is  low,  or  to 
an  orthogonal  isotropic  power  schedule  when  nominal  interference  is  high.  Furthermore,  the  transition  region  between  the 
latter  two  schedules  is  seen  to  be  very  sharp  in  terms  of  nominal  interference-to-noise  ratio  (INR).  For  MIMO  links  with 
correlated  channels,  the  corresponding  space-time  optimal  power  schedule  is  developed  based  on  the  knowledge  of  the 
channel  correlation  matrices.  It  is  shown  that  the  channel  correlation  has  a  great  impact  on  the  ergodic  capacity  and  the 
optimality  of  different  power  scheduling  approaches. 

14.  SUBJECT  TERMS 

N/A 

15.  NUMBER  OF  PAGES 

16 

16.  PRICE  CODE 

N/A 

17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION 

OR  REPORT  ON  THIS  PAGE  OF  ABSTRACT 

UNCLASSIFIED  UNCLASSIFIED  UNCLASSIFIED 

20.  LIMITATION  OF  ABSTRACT 

UU 

NSN  7540-0 1  -280-5500  Standard  Form  298  (Rev.2-89) 


Prescribed  by  ANSI  Std.  239-18 


298-102 


Enclosure  1 


686 


IEEE  TRANSACTIONS  ON  SIGNAL  PROCESSING,  VOL.  56,  NO.  2,  FEBRUARY  2008 


Space-Time  Power  Schedule  for  Distributed 
MIMO  Links  Without  Instantaneous  Channel  State 
Information  at  the  Transmitting  Nodes 

Yue  Rong,  Member,  IEEE ,  Yingbo  Hua,  Fellow,  IEEE ,  Ananthram  Swami,  Fellow,  IEEE ,  and 

A.  Lee  Swindlehurst,  Fellow,  IEEE 


Abstract — A  space-time  optimal  power  schedule  for  multiple 
distributed  multiple-input  multiple-output  (MIMO)  links  without 
the  knowledge  of  the  instantaneous  channel  state  information 
(CSI)  at  the  transmitting  nodes  is  proposed.  A  readily  computable 
expression  for  the  ergodic  sum  capacity  of  the  MIMO  links  is 
derived.  Based  on  this  expression,  which  is  a  non-convex  function 
of  power  allocation  vectors,  a  projected  gradient  algorithm  is 
developed  to  optimize  the  power  allocation.  For  a  symmetric  set 
of  MIMO  links  with  independent  identically  distributed  channels, 
it  is  observed  that  the  space-time  optimal  power  schedule  reduces 
to  a  uniform  isotropic  power  schedule  when  nominal  interference 
is  low,  or  to  an  orthogonal  isotropic  power  schedule  when  nominal 
interference  is  high.  Furthermore,  the  transition  region  between 
the  latter  two  schedules  is  seen  to  be  very  sharp  in  terms  of 
nominal  interference-to-noise  ratio  (INR).  For  MIMO  links  with 
correlated  channels,  the  corresponding  space-time  optimal  power 
schedule  is  developed  based  on  the  knowledge  of  the  channel 
correlation  matrices.  It  is  shown  that  the  channel  correlation 
has  a  great  impact  on  the  ergodic  capacity  and  the  optimality  of 
different  power  scheduling  approaches. 

Index  Terms — Multiple-input  multiple-output  (MIMO)  systems, 
space-time  power  schedule,  wireless  mesh  networks. 


I.  Introduction 

IN  a  large  wireless  mesh  network  of  many  multiple-input 
multiple-output  (MIMO)  nodes,  multiple  MIMO  links  must 
share  a  common  frequency  band  concurrently  to  ensure  a  high 
spectral  efficiency  of  the  whole  network  [  1  ]— [3] .  Developing  op- 
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timal  power  schedules  for  a  set  of  co-channel,  concurrent  and 
neighboring  MIMO  links  is  therefore  important. 

Power  scheduling  for  multiple  MIMO  links  has  been  studied 
in  [4]— [7] .  In  [4],  a  space-only  (i.e.,  time-invariant)  power 
schedule  is  presented,  and  an  iterative  algorithm  leading  to  the 
Nash  equilibrium  [8]  is  developed.  In  [5],  the  same  space-only 
criterion  is  used,  but  a  projected  gradient  algorithm  is  devel¬ 
oped  that  yields  a  better  result.  In  [6],  the  space-only  approach 
is  considered  without  channel  state  information  (CSI)  at  the 
transmitting  nodes.  In  [7],  a  space-time  power  schedule  is 
proposed  that  generalizes  the  approaches  used  in  [4]  and  [5], 

In  this  paper,  we  focus  on  networks  with  symmetric  topology, 
and  assume  that  the  network  topology  is  known  to  a  scheduler, 
and  each  receiving  node  knows  the  instantaneous  CSI  of  the 
link-of-interest.  For  a  given  set  of  transmitting-receiving  pairs, 
we  present  a  space-time  optimal  power  schedule  without  instan¬ 
taneous  CSI  at  the  transmitting  nodes.  This  work  goes  beyond 
that  in  [7],  which  assumes  instantaneous  CSI  at  the  transmit¬ 
ting  nodes  and  also  beyond  [6]  which  assumes  a  time-invariant 
transmitting  covariance  matrix  at  each  link. 

In  the  absence  of  instantaneous  CSI  at  the  transmitting  nodes, 
the  knowledge  of  some  statistical  properties  of  the  CSI  is  nec¬ 
essary  for  designing  optimal  power  schedules.  We  will  aim  to 
optimize  an  ergodic  network  capacity  with  respect  to  the  statis¬ 
tical  distributions  of  the  CSI.  In  other  words,  we  consider  a  “fast 
fading”  scenario  where  the  CSI  changes  rapidly  with  respect  to 
the  delay  requirement,  or  equivalently,  the  effective  data  rate 
that  we  are  interested  in  is  for  a  time  interval  during  which  the 
CSI  changes  many  times. 

We  will  consider  two  cases  of  the  statistical  properties  of  the 
CSI.  In  the  first,  the  channel  matrix  between  every  two  nodes 
consists  of  independent  and  identically  distributed  (i.i.d.)  com¬ 
plex  Gaussian  entries.  In  the  second,  each  channel  matrix  con¬ 
sists  of  correlated  entries.  The  first  case  will  be  referred  to  as  the 
case  of  i.i.d.  channels.  The  second  case  will  be  referred  to  as  the 
case  of  correlated  channels. 

An  important  result  in  this  paper  is  a  “closed  form”  expres¬ 
sion  for  the  ergodic  sum  capacity  of  multiple  MIMO  links.  This 
result,  following  from  the  quadratic  forms  of  Gaussian  random 
matrices  [9]  and  [10],  is  expressed  in  terms  of  finite  sums  and  a 
simple  one-dimensional  integral.  It  is  readily  computable.  An¬ 
other  important  result  in  this  paper  is  the  development  of  a  pro¬ 
jected  gradient  algorithm  that  allows  one  to  maximize  the  er¬ 
godic  sum  capacity  and  hence  to  compute  the  corresponding  op¬ 
timal  power  schedule. 
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To  illustrate  the  potential  of  the  space-time  optimal  power 
schedule,  we  will  compare  it  with  the  following  schemes. 

•  Uniform  isotropic  power  schedule'.  All  I  inks  conduct  trans¬ 
missions  at  the  same  time  and  the  same  frequency,  and  the 
source  covariance  matrix  used  by  all  links  is  proportional 
to  the  identity  matrix. 

•  Orthogonal  isotropic  power  schedule'.  All  links  conduct 
transmissions  in  orthogonal  channels  (such  as  different 
times  and/or  different  frequencies),  and  the  source  covari¬ 
ance  matrix  used  by  all  links  is  proportional  to  the  identity 
matrix. 

•  Low-rank  power  schedule'.  Each  link  uses  a  low  rank 
source  covariance  matrix. 

For  a  network  of  symmetric  MIMO  links  with  i.i.d.  chan¬ 
nels,  we  will  show  that  the  space-time  optimal  power  schedule 
consistently  outperforms  the  low  rank  power  schedule.  Inter¬ 
estingly,  the  space-time  optimal  power  schedule  reduces  to  a 
uniform  isotropic  power  schedule  when  nominal  interference 
level  is  low,  or  to  an  orthogonal  isotropic  power  schedule  when 
nominal  interference  level  is  high.  Furthermore,  the  transition 
between  the  latter  two  schedules  is  very  sharp  along  a  nom¬ 
inal  interference- to- noise  ratio  (I  NR)  axis.  All  space-only  power 
schedules,  including  those  in  [6],  are  shown  to  be  suboptimal 
compared  to  ours. 

It  is  important  to  note  that  by  I  NR,  we  refer  to  a  nominal 
INR  unless  specified  otherwise.  The  nominal  IN R  can  be  kept 
constant  while  the  actual  INR  changes,  which  will  be  further 
explained. 

We  will  also  show  that  the  channel  correlation  reduces  the 
network  capacity  when  nominal  interference  level  is  low,  but 
increases  the  network  capacity  when  nominal  interference  level 
is  high.  When  the  channel  correlation  is  known,  we  will  show 
that  the  eigenspace  of  the  optimal  source  covariance  matrices 
of  each  transmitting  node  is  equal  to  that  of  the  corresponding 
channel  covariance  matrix.  In  other  words,  the  channel  should 
be  “whitened”  by  the  source  covariance  matrices. 

The  rest  of  this  paper  is  organized  as  follows.  In  Section  II,  the 
system  model  of  distributed  MIMO  links  is  introduced,  and  the 
problem  formulation  is  given.  A  readily  computable  form  of  the 
ergodic  sum  capacity  of  the  M  IM  0  links  with  i.i.d.  channels  is 
derived  in  Section  III.  In  Section  IV,  we  present  the  space-time 
optimal  power  scheduling  scheme.  In  Section  V,  we  investi¬ 
gate  the  impact  of  the  channel  correlation  on  the  network  ca¬ 
pacity  and  study  the  corresponding  space-time  optimal  power 
schedule.  Numerical  examples  are  shown  in  Section  VI.  Con¬ 
clusions  are  drawn  in  Section  VII. 


II.  Problem  Formulation 

We  consider  data  transmissions  in  a  wireless  mesh  network 
of  M  IM  0  nodes.  We  focus  on  a  time  window  where  a  set  of  dis¬ 
tinct  pairs  of  transmitting  and  receiving  nodes  in  a  given  neigh¬ 
borhood  need  to  accomplish  their  data  transmissions.  Typically, 
a  pair  of  transmitting  and  receiving  nodes  consists  of  two  ad¬ 
jacent  nodes  for  power  and  bandwidth  efficiency.  We  also  con¬ 
sider  the  spectral  region  where  the  channel  frequency  response 
is  “flat”  (i.e.,  frequency- invariant).  The  key  problem  that  we  will 
address  in  this  paper  is  how  to  design  the  power  allocation  and 


the  source  covariance  matrix  for  each  transmitting  node  so  that 
the  network  throughput  during  this  time  window  is  maximized. 

There  are  two  simple  solutions.  One  is  to  force  all  transmit¬ 
ting  nodes  to  transmit  at  the  same  time  and  at  all  time,  and  the 
source  covariance  matrices  are  optimized  under  this  constraint. 
We  call  this  solution  the  space-only  power  schedule.  The  other 
solution  is  to  slice  the  time  into  multiple  slots,  and  during  each 
time  slot  only  one  transmitting  node  is  actually  transmitting. 
This  is  known  as  time-division  multiplexing  access  (TDM  A). 
With  TDM  A,  the  power  allocation  for  each  link  is  either  zero  or 
a  full  power  during  each  time  slot. 

There  is  a  more  general  solution  where  during  an  interval  of 
multiple  time  slots,  the  power  allocations  as  well  as  the  source 
covariance  matrices  for  all  links  are  treated  as  functions  of  the 
time  slots.  These  functions  can  be  optimized  to  maximize  the 
total  network  throughput.  This  is  the  approach  that  we  will 
follow  in  this  paper.  We  call  this  approach  space-time  optimal 
power  schedule.  Note  that  the  time  slots  here  can  be  replaced 
by  frequency  bins  or  a  combination  of  timeslots  and  frequency 
bins.  Both  time  and  frequency  are  temporal  characteristics  with 
respect  to  space.  But  for  a  simple  exposure  of  the  concept  of 
our  scheme,  we  will  only  refer  to  time  slots. 

For  convenience,  we  will  call  the  transmitting  node  and  the 
receiving  node  of  the  7th  link  as  the  7th  transmitting  node  and 
the  7th  receiving  node,  respectively.  Assuming  that  there  are  L 
links,  it  is  sufficient  to  design  the  optimal  power  schedule  over 
L  time  slots  [7],  We  will  assume  that  the  channel  matrices  re¬ 
main  constant  over  L  consecutive  time  slots  but  change  ran¬ 
domly  over  an  interval  of  many  multiples  of  L  time  slots.  We 
will  design  the  power  schedule  to  maximizean  ergodic  network 
capacity  which  is  averaged  over  the  statistical  distribution  of  the 
channel  matrices.  This  capacity  is  achievable  (approximately) 
over  the  interval  of  many  multiples  of  L  time  slots.  Therefore, 
ourtheory  is  valid  for  “fast  fading"  channels,  i.e.,  the  time  delay 
due  to  encoding  and  decoding  over  many  multiples  of  L  time 
slots  is  tolerable. 

We  assume  that  each  node  has  /V  antennas.  The  vector  of  the 
received  signal  at  the  7th  receiving  node  and  in  the  tth  time 
slot  can  be  written  as 

L 

y  i(t)  =  x/. ,  1 1;.  ,x ,  ( / )  +  V;(i). 

j=l,j^i 

(1) 

FI  ere,  denotes  the  average  gain  of  the  7th  transmitting  node 
and  the  7th  receiving  node,  which  is  for  a  link-of-interest. 
j  /  7,  is  the  average  gain  between  the  7th  transmitting  node 
and  the  7th  receiving  node,  which  is  for  an  interfering  link.  H;j 
is  the  instantaneous  TV  x  TV  channel  matrix  (normalized  by  the 
averaged  channel  gain)  between  the  jth  transmitting  node  and 
the  jr'th  receiving  node.  Xi(i),  7  =  denotes  the  TV  x  1 

vector  of  the  transmitted  signal  from  the  7th  transmitting  node 
within  the  slot  t.  Within  time  slot  t,  x,;(i)  is  a  random  function 
of  the  continuous  time  r  which  is  omitted.  v.;(i), 
is  the  TV  x  1  vector  of  the  i.i.d.  additive  (slot-invariant)  white 
Gaussian  noise  (AWGN)  with  zero  mean  and  unit  covariance 
matrix  CV!  =  Ijv-  A  nd  IiV  denotes  the  TV  x  TV  identity  matrix. 
Fig.  1  shows  an  example  of  a  circular  network  with  three  con¬ 
currently  communicating  M  IM  0  links.  A  list  summarizing  the 
notations  used  for  the  system  model  is  given  in  A  ppendix  A . 
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With  a  given  set  of  transmit  covariance  matrices  P;(f),  i,t  = 
L ... . .  /.,  the  averaged  ergodic  sum  capacity  of  L  links  over  L 
time  slots  can  be  written  as 


L  L 


Ia(P)  =  TEn  £El0& 


.  i=l  2=1 


L  N 


+/3i,lHM-Pz(f)H"zR-1(i) 


(2) 


where  |  •  |  denotes  the  determinant  of  a  matrix  or  the  absolute 
value  of  a  scalar,  £h{-}  stands  for  the  statistical  expectation 
with  respect  to  all  channel  matrices  H  =  [H^1? . . . 
and  (-)T  denotes  the  matrix  transpose.  In  (2),  P  is  a  matrix 
stacking  the  source  covariance  matrices  of  all  links: 

p  =  [PJ,...,ptl]t 

Pi^  [pf(i),...,p J(L)]T  /-I...,/. 

and  R,  ( t )  i  s  the  (actual )  i  nterference- pi  us- noi  se  covari  ance  ma¬ 
trix  at  the  2th  receiving  node  in  the  7th  time  slot 


Fig.  1.  Circular  wireless  network.  Solid  circles:  source  nodes:  hollow  circles: 
destination  nodes.  Solid  lines:  signal  streams:  dashed  lines:  interference  streams. 


Thefirstterm  in  (1)  isthe  signal  of  interest  at  the  2th  receiving 
node,  and  the  second  term  isthe  sum  of  interfering  signalsfrom 
all  otherW- 1  transmitting  nodes.  We  assume  that  all  the  trans¬ 
mitted  signals  are  Gaussian  distributed  with  zero  mean  and  their 
covariance  matrix  is  periodically  time  varying  with  the  period 
equal  to  L  timeslots.  We  will  write  Pz(7)  =  E  {xi(f)xi(i)H}, 

/  •  -.i _ ,L,  where  E{-}  stands  for  statistical  expectation,  and 

(■)H  denotes  matrix  Hermitian  transpose.  M  oreover,  we  assume 
that  (l/L)  Y^t=i  MPiW}  =  Pi,  i  =  where  P{  is 

the  averaged  transmit  power  from  the  2th  transmitting  node,  and 
tr{-}  stands  for  the  trace  of  a  matrix.  We  will  need  further  as¬ 
sumptions  as  shown  below. 

•  There  is  no  coding  cooperation  between  interfering  trans¬ 
mitting  nodes.  Namely,  the  source  information  is  encoded 
independently  at  different  transmitting  nodes.  Cooperative 
coding  among  interfering  M  IM  0  links  would  significantly 
increase  the  networking  complexity,  which  is  undesirable. 

•  The  entries  of  Hzj,  2,  j  =  1, ...  ,L,  are  complex  Gaussian 
distributed  with  zero  mean  and  unit  variance.  (As  men¬ 
tioned  before,  H iv7-  remains  constant  over  L  time  slots  but 
varies  randomly  during  an  interval  of  many  multiples  of  L 
time  slots.  For  a  stationary  network,  such  random  fading 
can  be  induced  by  randomly  changing  the  phase  of  each 
transmitting  antenna.)  Thus,  the  averaged  signal-to-noise 
ratio  (SNR)  at  each  antenna  of  the  2th  receiving  node  is 
given  by  Ai  =  / %iPi ,  while  the  nominal  INR  from  the  jth 
transmitting  node  to  each  antenna  of  the 2th  receiving  node 
is  defined  as  B.hl  =  iijjPj.  We  treat  pitj,  2  =  1 , ,22 
as  deterministic  parameters  both  in  the  analysis  and  the 
simulations. 

•  There  is  no  instantaneous  CSI  knowledge  at  any  transmit¬ 
ting  node,  and  the  2th  receiving  node  only  knows  the  CSI 
of  the  link-of-interest,  i.e.,  HZii. 


L 

R i(t)=  Y,  ■^H,(P,(/!H,/.//  •  Iv. 
i  U  :i/:  i 

in  this  paper,  we  aim  at  maximizing  the  ergodic  capacity 
(2)  by  choosing  the  transmit  covariance  matrices  P;(f),  2,7  = 
1  This  problem  can  be  compactly  formulated  as 

max  Ia  (P)  (3) 

p 

1  L 

s.t.  -£tr{P  i(t)}  =  Pi  (4) 

^  t= 1 

Pi(t)  >  0  i,t  =  i,...,L  (5) 

where  (4)  represents  the  average  transmit  power  constraint  at 
each  transmitting  node.  Note  that  in  this  paper,  we  focus  on  net¬ 
works  with  a  symmetric  topology.  Therefore,  no  link  suffers  a 
fairness  problem  under  the  ergodic  sum  capacity.  W  hen  the  net¬ 
work  is  asymmetric,  fairness  issues  arise.  Then,  the  cost  func¬ 
tions  must  be  changed.  One  commonly  used  criterion  is  known 
as  “proportional  fairness,”  which  can  be  formulated  as  sum  of 
logarithmsof  link  capacities.  The  gradient-based  method  can  be 
applied  similarly.  But  the  details  would  go  beyond  the  scope  of 
this  paper. 


III.  Ergodic  Sum  Capacity  for  MIMO  Links 
W ith  i.i.d.  Channels 

in  this  section,  we  derive  a  closed-form  expression  for  the  er¬ 
godic  sum  capacity  (2)  for  the  case  of  i.i.d.  channels.  The  corre¬ 
lated  channel  casewill  be  addressed  in  Section  V .  Forsimplicity, 
we  first  derive  a  closed-form  expression  for  the  following  er¬ 
godic  sum  capacity  at  one  time  slot: 


/(p1(i),.^.,pL(i))  =  ^i£;iog2 


.  i— 1 


LAT 


-W;,;H;.;P;(7)H''Rrl(/) 


(6) 
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We  will  later  generalize  the  result  to  obtain  a  closed-form  ex-  where 
pression  for  the  ergodic  sum  capacity  averaged  over  all  time 
slots  as  shown  in  (2). 

Let  us  denote 


P  i(t)  =  f(t)  (7) 

as  the  eigendecomposition  of  P i(t),  where  Ui(i)  is  the  unitary 
matrix  of  eigenvectors,  and 


*i(t)=  E  / (Z  )H •  j  -  T  v 

^  =  [Hi; Hi,i] 

=  [H^i, .  .  .  ,  1, 

A i(t)  =  diag  | 

A i(t)  =diag{ 


■  ■  ■  5  Ai;i_i(i), 


D  i(t)  =  diag{dii(f),  di2(t), ....  diN(t)} 

is  an  NxN  diagonal  matrix  of  all  eigenvalues.  For  convenience, 
we  will  use  the  N  x  1  column  vectors  di(i),  i  =  1 
defined  as 

d  i(t)  =  [dn(t),  di2(t ), . . . ,  diN(t)]T. 


Kj(t)  =Pi,jdj(t). 

It  can  be  seen  from  (10)  that  the  ergodic  sum  capacity  ex¬ 
pression  is  a  summation  of  2 L  logarithmic  terms  all  having  a 
similar  structure.  For  convenience,  we  can  now  focus  on  one  of 
the  terms 


Substituting  (7)  into  (6),  we  can  rewrite  (6)  as 


m 


. .. ,  dL(t))  =  \  E  lo§2 

U=i  I 
+/?i,.iHi,i(f)Dl(i)Hfi(f)Rr1(i) 


where 


(8) 


H(*)  = 


=Hi,iUi(i) 
H  =  HijUj(f) 


L 

R 4(t)=  E  ■KiU,iu)D,u)ii!'iii)  ■  iv. 

SinceHij,7,j  =  1  has  i.i.d.  Gaussian  entries  and  Ui(i) 
is  unitary,  the  statistics  of  H are  identical  to  those  of  H;j 
[11],  H ereafter,  for  notational  simplicity,  we  drop  the  tildes  on 
the  H  matri  ces  and  the  dependence  of  H  matri  ces  on  t,  and  w  rite 
the  ergodic  sum  capacity  expression  (8)  as 


7(d1(t),...,dL(t)) 


=  77h  E  l0S2  l1*  +  /5t,iHi,iDi(t)H^RT1(i)| 
U=i 

=  E  {log2  |R i(t)  +  A,iHi,iDi(i)Hfi|} 


i=l 


EH{log2  |Ri(i)|}J 


=  E  (  (l0S2  |liV  +  H,:A.,(i)Hf  |  } 


i= 1 


-^{logzIlrv  +  HA^fifl} 


(9) 


(10) 


iii(Ai(i))  =  Eh,  {log2  \IN  +  H,A.,(*)Hf|}.  (11) 

It  can  be  seen  that  (11)  is  equivalent  to  the  ergodic  capacity  of  a 
point-to-point  MIMO  system  under  AWGN  with  NL  transmit 
antennas  and  N  receive  antennas.  A  closed-form  expression 
for  the  ergodic  capacity  of  such  a  system  has  been  shown  in 

[9] ,  where  a  determinant  representation  for  the  distribution  of 
quadratic  forms  of  a  complex  Gaussian  matrix  [10]  has  been 
used.  Using  the  results  in  [9],  (11)  can  be  written  as 

N—l  NL 
=  log2(e)  EE 

n= 0  fc=l 

where  [see  (12)— (14),  shown  at  the  bottom  of  the 
next  page],  (Note  that  there  are  typos  in  (16)  of  [9], 
where  the  terms  in  the  first  summation  should  be 
(z!(-l)(i^ty(*  -  r)!)7(:fifl^e1/75,i  (1/7),  using  the  same 
notations  as  in  [9]  except  for  Si.)  H ere,  =  [A^(i)]fc;fc, 
k  =  1,...,NL,  denotes  the  (k,k) th  element  of  matrix  A, (t), 
and  Si  (a;)  =  /“  e~T /rdr  is  the  exponential  integral  function 
of  order  1  [12,  eq.  5.1.1].  It  is  worth  noting  that  (12)— (14)  are 

derived  under  the  assumption  that  all  Xt^k,  k  =  I . N /,, 

have  distinct  values.  W  hen  some  of  them  are  identical,  we  can 
obtain  the  closed-form  ergodic  sum  capacity  expression  by 
deriving  the  limit  of  (11)  with  respect  to  those  common  values 
of  \i,k  using  L'H ospital’s  rule.  For  numerical  evaluation,  it 
is  sufficient  to  slightly  and  randomly  perturb  these  identical 
values  of  Xt^k  because  all  functions  are  continuous  and  Xt^k 
are  deterministic. 

Similar  to  (11),  we  can  write  each  of  the  2 L  summands  in 

(10)  in  a  closed  form.  Thus,  the  ergodic  sum  capacity  can  be 
represented  as 


/(di(f),..t,dL(f)) 


L 

=  log2(e)^ 
i= 1 


~N- 1  NL 

EE  CtiknQijI'i  ^t,i,k) 

_n=0  k= 1 
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N-lN(L-l) 

-E  E  'LHiknQ  (tL*  \t,i,k) 

n=0  k=l 


(15) 


where  =  [A;(f)j  kk,  A-  1),  ctifen  is  given 

by  (12)  and  (13),  and  uukn  is  given  by  (16)  and  (17),  shown  at 
the  bottom  of  the  page. 

The  result  (15)  can  be  readily  used  to  write  the  averaged  er- 
godic  sum  capacity  (2)  in  closed  form.  First,  we  rewrite  (2)  as 
(18)  and  (19),  shown  at  the  bottom  of  the  page.  Then,  by  fol¬ 
lowing  the  steps  given  above,  we  can  write  (19)  in  closed  form 
as 


T>(d) 


log2(e) 


EE 

t  i  i  i 


TV— 1  NL 


EE  C-tiknQijl')  \ 


1 t.i.k  ) 


,n= 0  k= 1 


N-lN(L-l) 

-E  E  'U'tiknQijI'i  X t,i,k ) 

?i.=0  k= 1 


(20) 


where  d  is  a  vector  stacking  the  power  scheduling  parameters 
of  all  links  at  all  L  time  slots 


d;  =  [df(l),...,df  (L)]t 


7  =  1 ...  ,L 


and  dj  is  a  vector  stacking  the  power  scheduling  parameters  of 
the  7th  link.  In  (20),  ctikn  and  uukn  can  be  written  using  (12) 
and  (13),  and  (16)  and  (17),  respectively. 

As  shown  in  (20),  theergodic  sum  capacity  is  now  expressed 
as  a  finite  summation  involving  rational  functions  and  exponen¬ 
tial  integration  functions  of  the  power  scheduling  vectors  d; (A), 
i,t  =  1, . .. ,L,  of  all  transmitting  nodes.  The  exponential  inte¬ 
gration  function  is  available  in  many  software  such  as  M  ATLAB 
and  M  athematica.  Thus,  (20)  is  easy  to  compute.  For  a  better  un¬ 
derstanding,  the  function  Q(n,  A)  is  plotted  in  Fig.  2  versus  A 
for  n  =  0,  1,  2,  3.  The  closed-form  expression  (20)  enables  us 
to  numerically  optimize  the  averaged  ergodic  sum  capacity  and 
hence  the  power  scheduling. 

The  accuracy  of  the  ergodic  sum  capacity  expression  (20)  is 
verified  by  Fig.  3.  This  figure  shows  both  the  analytical  and 
the  simulation  results  of  the  ergodic  capacity  versus  the  nom¬ 
inal  INR  (i.e.,  B)  of  a  symmetric  network  (B^j  =  B),  with 
L  =  2  active  links,  where  each  node  has  TV  =  2  antennas.  FI  ere, 
we  set  Pi  =  20  dB ,  0iti  =  1,7  =  1, . . . ,  L,  and  0ij  =  0, 
i,j  =  1, . . .  ,77,  j  ^  7.  The  ergodic  capacity  shown  is  a  per-link 
ergodic  capacity  which  is  the  ergodic  sum  capacity  divided  by 
the  number  of  active  links  L.  For  the  simulation  results,  1000 
channel  realizations  were  used.  A  randomly  selected  power  al¬ 
location  vector  was  used  for  each  link.  It  can  be  seen  from  Fig.  3 
that  the  analytical  result  is  consistent  with  the  simulation  result. 


Ctikn  — 


n! 


vh^k 


Vtikn 


btikn  — 


jr^k 

Xy  •  •  ♦  ^t,i,jN-n- 1  5  Ti  =  0,  .  .  .  ,  N  —  2 

i<ji<-<jN-n-i<NL 

1,  n  =  N  —  1. 


poo 

Q(n,Xt,i,k)=  /  ln(l  +  x)xne~^x/Xi’i’k'ldx 

(4)+eee  "!(-1)<”',a«2 

r=0  T>  ■  \At,z,kJ  r=l  s=0  h= O 


(n  —  r)\(r  —  s  —  1  —  h)\(r  —  s) 


(12) 

(13) 


(14) 


'U'tikn  —  ' 


bfikn  — 


^  \  ^  N— n— 1  yN(L— 2) — 1 


n! 


jTJ  ^ 

vh^k 


-1 


t.i.h  I  I  Otikn 


jr^k 

^ t,i,j i  •  •  •  Xt,i,j N,_in_1  ■  n  =  O,. . .  ,N  —  2 

l<Jl<"'<JJV-n-l<JV(L  —  1) 

1,  n  =  N  —  1. 


(16) 

(17) 


L  L 


MP)=tEE  log2 


t=l  2=1 
L  L 


Ri«)  + 


-77H{log2  |Ri(i)|} 


=  7EEPh,  {log2  |Ijv  +  H,Ai(A)Hf  |}  -  {log2  |lv  +  H;A;(A)Hf  |}  . 


(18) 


(19) 


t=l  2=1 
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Fig.  2.  Function  Q(n,  A)  for  n  =  0, 1,  2,  3. 


where  (22)  is  the  set  of  the  transmit  power  constraints  at  all 
transmitting  nodes,  and  ||  -  Hx  denotes  the  sum  norm  (or^  norm) 
of  a  vector.  For  a  vector  x,  x  >  0  means  that  each  entry  of  x  is 
nonnegative. 

The  results  in  [6]  show  that  when  the  IN  R  is  sufficiently  low, 
theergodic  sum  capacity  (15)  isaconcavefunctionof  the  power 
allocation  vectors  di(i),. . .  ,dL(i),  but  when  the  IN R  is  suffi¬ 
ciently  high,  (15)  becomes  a  convex  function  of  the  power  allo¬ 
cation  vectors.  (Si nee  the  power  allocation  in  [6]  is  independent 
of  time,  the  IN R  in  [6]  is  an  instantaneous  INR,  an  averaged 
INR  and  also  the  nominal  INR.)  However,  in  general,  itcan  be 
seen  from  (10)  that  due  to  the  mutual  interference  among  dif¬ 
ferent  links,  the  ergodic  sum  capacity  is  neither  a  convex  func¬ 
tion,  nor  a  concave  function,  of  the  power  allocation  vectors 
d !(i), . . . ,  dL(t).  Similarly,  (20)  is  neither  a  convex  nor  a  con¬ 
cave  function  of  the  power  scheduling  vector  d  in  the  general 
INR  region.  Thus,  in  general,  (21)— (22)  is  a  nonconvex  opti¬ 
mization  problem. 

The  objective  function  (21)  is  quite  complicated  and  does  not 
have  a  clear  structure.  The  gradient-based  methods  are  effec¬ 
tive  for  such  a  problem  [13,  ch.  1],  These  methods  guarantee 
the  convergence  to  a  stationary  point,  provided  that  proper  step 
sizes  are  chosen.  M  oreover,  si  nee  the  constraints  (22)  are  simple 
linear  constraints,  the  projected  gradienttechnique  [13,  sec.  2.3] 
can  be  applied  to  solve  the  problem  (21)— (22).  We  first  give  a 
brief  introduction  to  the  projected  gradient  algorithm. 

Let  us  consider  the  following  constrained  optimization 
problem 

max/(x)  s.t.  xeT 

where  /(•)  is  a  continuously  differentiable  scalar  function,  and 
A"  is  a  nonempty,  closed,  and  convex  set.  The  projected  gradient 
algorithm  is  a  feasible  direction  method.  The  algorithm  starts  at 
an  initial  point  #.  At  the  fcth  iteration,  x^fc)  is  updated  as 


Fig.  3.  Comparison  of  the  analytical  and  simulation  results  of  the  per-link  er¬ 
godic  capacity  versus  nominal  INR  for  the  case  of  i.i.d.  channels.  Here,  N  =  2, 
L  =  2. 


IV.  Space-Time  Power  Schedule 

It  has  been  shown  in  [7]  that  by  applying  a  space-time  power 
schedule  where  the  source  covariance  matrices  are  allowed  to 
be  functions  of  time,  a  larger  (averaged)  ergodic  sum  capacity 
can  be  achieved.  However,  the  power  scheduling  scheme  in  [7] 
requires  the  instantaneous  CSI  knowledge  at  the  transmitting 
nodes,  which  is  valid  for  slow  fading  channels.  In  this  section, 
we  optimize  the  space-time  power  schedule  for  fast  fading 
channels  by  maximizing  the  ergodic  sum  capacity  expression 
derived  in  Section  III  for  which  the  instantaneous  CSI  is  not 
required  at  the  transmitting  nodes. 

Taking  the  power  constraint  of  each  active  link  into  account, 
our  space-time  power  scheduling  approach  becomes  the  fol¬ 
lowing  optimization  problem: 

max  Ia  (d)  (21) 

d 

s.t.  1 1 dy 1 1 x  =  LPi,  dj  >  0.  (22) 


x(fc+i)  =  x(fe)  +  6k  ^x(fc)  _  x(fc)  j  (23) 

where  Sk  e  (0, 1]  is  a  step  size,  and 


x(fc) 


=  proj 


sfcV/  (x«) 


(24) 


Here,  proj[-j  denotes  the  projection  onto  the  feasible  set  X,  and 
sk  is  a  positive  scalar.  From  (23)  and  (24),  we  can  see  that  to 
obtain  the  vector  x(fc),  we  first  take  a  step  sfcV /  (x(fc))  along 
the  gradient,  similar  to  the  steepest  descent  algorithm.  Then  we 
project  the  result  x(fc)  +  sfcV /  (#)  onto  the  feasible  set  X. 
It  has  been  shown  that  after  a  sufficient  number  of  iterations, 
the  sequence  x<fc)  converges  to  a  local  optimal  point  Xopt  [13, 
sec.  2.3], 

Two  major  issues  of  a  projected  gradient  algorithm  are:  1) 
computing  the  gradient  of  the  objective  function,  and  2)  pro¬ 
jecting  the  gradient-modified  vector  into  the  feasible  set. 

The  derivation  of  the  gradient  of  (20)  is  tedious  but  straight¬ 
forward.  The  complete  expression  for  the  gradient  of  (20)  with 
respect  to  dlm{t),l  =  1, . . .  ,L,  m  -  1, . . .  ,7V,  t  =  ij*. .  ,L  is 
listed  in  Appendix  B. 
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In  the  sequel,  we  show  how  to  carry  out  the  projection  op¬ 
eration.  Let  us  take  the  7th  link  as  an  example,  and  denote  the 
NL  x  1  column  vector  of  gradient  as 

A  f  a/g(d)  a/n(d)  lT 
St  lddil(i),'",ddiN(L)_  ' 

Then  taking  a  step  along  the  positive  gradient  (24),  the  ALL  x  1 
power  vector  is  updated  as 

dj  =  di  +  sgi. 

The  next  step  of  the  projected  gradient  algorithm  is  to  project 
d;  onto  the  feasible  region  of  power  vectors  defined  by  (22). 
In  fact,  the  projection  operation  can  be  seen  as  searching  for  a 
point  dj  in  the  region  of  (22),  which  has  a  minimal  Euclidean 
distance  to  the  point  d;.  Thus  for  the  7th  link,  the  projection 
operation  can  be  formulated  as  the  following  constrained  opti¬ 
mization  problem: 

A  2 

min  di  —  dj  (25) 

d, 

s.t.  d i  =  LPi ,  di  >  0,  7  =  1, . . . ,  L  (26) 

l 

where  ||  •  ||  denotes  the  Euclidean  norm  of  a  vector. 

The  problem  (25)— (26)  can  be  efficiently  solved  by  using  the 
Lagrange  multiplier  method  [14,  sec.  5.1].  Each  element  of  d; 
can  be  written  as 

din=[din-^\+ :  n  =  1, . . .  ,NL  (27) 

where  //  is  the  Lagrange  multiplier,  and  for  a  real  scalar  x, 

[x]  r  =  inax{./:.  ()}. 

The  Lagrange  multiplier  n  can  be  obtained  by  substituting 
(27)  back  into  (26).  We  have  the  following  equation: 

NL 

y:  max  { din  _  f  ’  °}  =  Lpi-  <28) 

n= 1  ~ 

The  left-hand  side  of  (28)  is  a  piecewise  linear  function  and 
monotonically  decreasing  with  respect  to  yt.  Thus,  (28)  can 
be  easily  solved  by,  for  example,  the  bisection  method  [14, 
pp.  145], 

N  ote  that  under  the  space— ti  me  power  schedul  i  ng  scheme,  the 
iterative  computation  of  the  power  allocation  vectors  for  all  links 
can  be  done  concurrently  or  sequentially.  The  concurrent  iter¬ 
ations  must  be  done  at  a  single  place  (centralized).  But  the  se¬ 
quential  iterations  can  be  done  in  a  distributed  manner,  i.e.,  each 
link  updates  its  own  power  allocation  vector  at  each  iteration.  At 
the  end  of  each  sequential  iteration,  all  I  inks  must  exchange  their 
updated  power  allocation  vectors  before  the  next  sequential  it¬ 
eration  takes  place.  Due  to  space  limitation,  weonly  present  the 
centralized  technique  in  detail  in  this  paper. 

At  the  fcth  iteration  of  the  centralized  technique,  the  power 
scheduling  vector  is  updated  as 

d./fc+1)  =  d/fe)  +  Sk  (d/fe)  -  d/fe))  (29) 


TABLE  I 

Procedure  of  the  Projected  Gradient  Power 
Allocation  Approach 

1)  Initialize  power  allocation  vectors  d,. 

Choose  step  sizes.  Set  k  =  0. 

2)  Set  k  :=  k  +  1. 

Calculate  the  gradient  of  (20)  using  (53)-(59). 

Let  d =  di(fc)  +skg(ik). 

Project  to  obtain  d)fc)  using  (27). 

Update  d,',;>  using  (29)-(30). 

3)  If  convergent,  end. 

Else  go  to  step  2. 

d;(^  =  proj  d;(fc)  +  skg\k)  ,  i  =  l,...,L  (30) 

where  Sk  and  sk  are  scalars  of  step  size.  They  can  be  chosen, 
for  example,  according  to  the  Armijo  rule  along  the  feasible  di¬ 
rection  [13,  pp.  225-226],  In  this  rule,  sk  =  s  is  a  constant 
throughout  the  iterations,  and  Sk  =  6mk,  where  mk  is  the  min¬ 
imal  nonnegative  integer  that  satisfies  the  following  inequality: 

r„  (d<*+1>)  -  h  (aw) 

>^~E(gP)T(di(‘,-d,w). 

2=1 

Here,  a  and  6  are  constants.  According  to  [13,  pp.  225-226], 
usually  a  is  chosen  close  to  0,  and  a  proper  choice  of  6  is  from 
0.1  to  0.5. 

The  steps  of  (29)  and  (30)  are  carried  out  for  all  L  links  and 
continue  till  the  sequence  d(fc)  converges.  Such  convergence  is 
guaranteed  according  to  the  principle  of  the  projected  gradient 
algorithm  using  the  A  rmijo  rulealong  thefeasible  direction  [13, 
Proposition  2.3.1],  The  convergence  criterion  is  taken  as 

maxabs  |d^’+1^  —  d^|  <  e 

where  maxabs{-}  denotes  the  maximal  absolute  value  among 
all  elements  of  a  vector,  and  e  is  a  positive  constant  close  to  0. 
Since  the  projected  gradient  algorithm  only  converges  to  a  local 
optimal  solution,  multiple  initializations  can  be  used  to  improve 
the  performance.  In  Section  VI,  we  will  study  the  achievable 
ergodic  capacities  with  different  number  of  initializations. 

The  procedure  of  applying  the  projected  gradient  technique 
to  solve  the  problem  (21)— (22)  is  summarized  in  Table  I.  Note 
that  there  may  be  peak  transmit  power  constraint  in  some  prac¬ 
tical  systems  such  that  di(t)  <  <7max,  where  dmax  denotes  the 
peak  transmit  power.  The  proposed  projected  gradient  approach 
can  be  easily  extended  to  solve  the  problem  (21)— (22)  with  the 
additional  peak  transmit  power  constraint. 

The  proposed  space-time  power  schedule  requires  only  the 
knowledge  of  averaged  SNR  and  nominal  INR  of  each  link  at 
the  transmitting  nodes.  This  knowledge  can  be  easily  obtained 
by  exploiting  the  topology  of  the  wireless  networks  and  the  av¬ 
eraged  transmit  power  of  each  transmitting  node,  which  in  turn, 
can  be  obtained  by  having  each  node  report  its  current  location 
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and  power  budget  to  the  scheduler  periodically.  I  n  practice,  the 
optimization  procedurecan  be  run  offli  ne  for  different  combina¬ 
tions  of  averaged  SN  R  and  nominal  IN  R.  The  resulting  optimal 
parameters  can  be  tabulated.  Then,  in  real-time  applications,  we 
only  need  to  look  up  this  table  to  select  the  optimal  power  pa¬ 
rameters.  It  is  worth  noting  that  although  the  instantaneous  SNR 
and  the  instantaneous  IN  R  may  vary  from  time  slot  to  timeslot, 
the  averaged  SN  R  and  the  nominal  IN  R  are  static  for  many  mul¬ 
tiples  of  L  timeslots  according  to  the  assumption  in  Section  II. 

Before  finishing  this  section,  we  want  to  discuss  a  special 
case  of  the  proposed  space-time  power  schedule.  When  only 
one  time  slot  is  considered  for  power  scheduling,  we  have  the 
space-only  power  schedule,  which  can  be  written  as  the  fol¬ 
lowing  constrained  optimization  problem: 

max  /( di,...,dL)  (31) 

s.t.  I  d;|  I  -  P,.  d;  >  0,  i  =  l,...,L.  (32) 

it  has  been  shown  in  [6]  that  at  a  sufficiently  low  interference 
level,  the  space-only  optimal  power  schedule  is  a  uniform 
isotropic  power  schedule  where  all  links  use  the  same  source 
covariance  matrix  and  the  source  covariance  is  the  identity  ma¬ 
trix.  While  at  a  sufficiently  high  interference  level,  it  is  shown  in 
[6]  that  the  space-only  optimal  power  schedule  becomes  a  low 
rank  power  schedule  where  each  link  uses  a  low  rank  source 
covariance  matrix.  However,  the  work  [6]  does  not  provide  a 
good  answer  for  the  intermediate  region  of  interference.  Our 
optimization  based  on  a  single  time  slot  yields  the  space-only 
optimal  power  schedulefor  any  given  interference  level,  which 
will  be  shown  in  Section  VI. 

Also  note  that  for  the  space-only  power  schedule,  the  nom¬ 
inal  interference  referred  to  in  this  paper  is  also  the  actual 
interference. 


V.  Channel  Correlation 

I  n  Sections  1 1 1  and  I V ,  we  assumed  that  the  channel  matrix  has 
i.i.d.  Gaussian  elements.  Yet,  in  some  practical  applications,  the 
elements  of  a  channel  matrix  can  be  correlated  due  to  reasons 
such  as  the  short  spacing  between  antenna  elements,  and/or  lack 
of  scatters  around  the  antennas  [15],  [16],  In  this  section,  we 
investigate  the  impact  of  channel  correlation  on  theergodic  sum 
capacity,  and  also  study  the  corresponding  space-time  optimal 
power  schedule.  A  s  we  aim  to  study  the  optimal  power  schedule, 
the  channel  correlati  on  at  the  transmitter  si  de  i  s  of  great  i  nterest. 
A  widely  accepted  model  for  such  a  channel  correlation  is  given 
by  [15],  [16] 


H-,j  =  H";,  O'/2  (33) 

where  Hf  •  has  i.i.d.  complex  Gaussian  distributed  entries  with 

',')j  ' 

zero  mean  and  unit  variance,  and  >  0  characterizes  the 
channel  correlation  at  the  transmit  side.  Although  the  M  IM  0 
channel  realizations  H,  ,  can  be  time  varying  from  slot  to  slot, 
we  assume  that  the  channel  statistics  0,  is  static  for  many  mul¬ 
tiples  of  L  time  slots. 


Let  us  denote 


A  v,r,v"  (34) 

astheeigendecomposition  of  the  channel  covariance  matrix  fi,, 
where  Vj  is  a  unitary  matrix  of  eigenvectors,  and  Tj  is  a  di¬ 
agonal  matrix  containing  all  the  eigenvalues.  Substituting  (34) 
back  into  (33),  we  have 

H  —  H"V  r  ;>v,/ 

—  ',:i  v  JL  j  v  j  • 

Let  us  introduce  the  following  new  matrices 

S;., A  n,,v(rjl/2  -  n;;;,v,  (35) 

Qj  (*)  =  r/2  vf  p;/  (f )  V,T'j/2 .  (36) 

it  can  be  seen  from  (35)  that  since  V7  is  a  unitary  matrix,  S itj 
has  the  same  di  stri  buti  on  as  H^- [  1 1] .  I  n  other  w ords,  the  entri  es 
of  S t.j  are  i.i.d.  Gaussian  distributed  with  zero  mean  and  unit 
variance.  Substituting  (35)  and  (36)  back  into  (6),  we  have 

I(Qi(f),  •  •  • ,  Q L(t))  =  ES  {  £  log2  IN 

l  i= j, 

+/3i,LSi,iQl(i)SflT-1(i)  |  (37) 
where  S  =  [S^,...,S^JT,  and 


T i(t)=  ]T  7,,S;.(Q;(/)S"  •  Iv. 

When  the  channel  covariance  matrix  0;  is  unknown  at  the 
transmitting  nodes,  a  straightforward  (suboptimal)  space-only 
power  schedule  is  P i(t)  =  Pi/NIN.  By  substituting  (36)  into 
(37),  we  have 

/(lr»-""xrA')=£s(Y;io&iiv 

+^sl,;r;sflF-1 1  (38) 

where 


I’’;-  E  ^SvTjSZ  +  Is-  (39) 

From  (38)  and  (39),  we  can  see  that  in  this  case,  the  ergodic 
sum  capacity  is  equivalent  to  that  of  the  multiple  MIMO  links 
with  i.i.d.  fading  channel  and  the  source  covariance  matrices 
Pi/NTi,  i  =  L.  Thus,  the  impact  of  the  channel  cor¬ 
relation  on  theergodic  sum  capacity  for  the  space-only  case  can 
be  analyzed  as  follows. 
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Fig.  4.  Comparison  of  ergodic  capacities  of  M IM  0  networks  with  channel  cor¬ 
relation  a)  =  (0,0),  u  =  (0.4, 0.5 ),  u,-  =  (0.8,0.9),  ora;  =  (1,1). Here, 
N  =  2,  L  =  2. 


When  the  interference  level  is  low  such  that  it  can  be  ne¬ 
glected,  we  have  L  parallel  point-to-point  MIMO  links.  It  is 
well  known  [11]  that  for  a  single  link  with  unknown  CSI  at 
the  transmitter,  the  optimal  source  covariance  matrix  is  a  scaled 
identity  matrix.  If  the  channel  is  correlated,  P;  /  clN,  where  c 
is  a  scalar,  then  the  ergodic  capacity  is  reduced  by  using  as 
the  source  covariance  matrix.  For  the  case  of  highly  correlated 
channels,  the  values  of  the  diagonal  elements  of  V;  are  widely 
spread.  This  may  lead  to  a  severe  performance  degradation  of 
the  ergodic  sum  capacity. 

When  the  interference  level  is  high,  it  has  been  shown  in  [6] 
that  the  optimal  space-only  power  schedule  is  equivalent  to  a 
low  rank  power  schedule.  W  hen  thechannel  is  highly  correlated, 
r;  has  the  tendency  to  be  rank-deficient.  In  this  case,  thechannel 
correlation  “assists" alow  rank  power scheduleand  hencecould 
improve  the  ergodic  sum  capacity. 

Fig.  4  shows  a  numerical  example  of  the  impact  of  channel 
correlation  on  the  per-link  ergodic  capacity  of  distributed 
MIMO  links  with  L  =  2  active  and  symmetric  links,  where 
each  node  has  IV  =  2  antennas  and  the  space-only  optimal 
power  schedule  is  applied.  FI  ere,  we  set  p  =  20  dB,  0iti  =  1, 
*  =  and  0id  =  0,  i,j  =  1  The 

exponential  correlation  model  has  been  adopted  to  simulate  the 
channel  correlation,  i.e.,  the  (fc,Z)th  element  of  matrix  Q;  is 
given  by 


[ft*ki=4fc  l\  0<wi<l, 

i  =  l7...,L,  k,l  —  1 _ ;Y.  (40) 

We  will  use  w  =  (04,., ..  ,ljl)  to  capture  the  channel  corre¬ 
lation  of  all  L  links.  Fig.  4  supports  our  analysis,  i.e.,  the  er¬ 
godic  capacity  decreases  as  the  channel  correlation  increases 
when  the  interference  is  low,  but  the  ergodic  capacity  increases 
as  the  channel  correlation  increases  when  the  interference  is 
high.  In  particular,  the  ergodic  capacity  becomes  very  robust 
against  I  NR  when  the  channels  are  fully  correlated,  i.e.,  when 
w  =  (1,1). 


In  the  sequel  of  this  section,  we  develop  the  space-time  op¬ 
timal  power  schedulefor  MIMO  links  with  correlated  channels. 
We  assume  that  all  channel  covariance  matrices  of  all  links  are 
known  at  each  transmitting  node. 

The  overall  network  ergodic  sum  capacity  averaged  over  L 
time  slots  can  be  written  as 

mq)=HeeioS2 

U=i  i= 1 

where  Q  is  a  matrix  stacking  the  Q  matrices  of  all  links: 

Q  =  [Qf,...,Ql]T 

Q;  =  [Qf  (1),  (L)] T 

The  space-time  optimal  power  schedule  can  beformulated  as 
the  following  constrained  optimization  problem: 

max  7a(Q) 

Q 

s.t.  ^Etr{rr1/2 Qi(*)rr1/2}  =  ^ 

n  t= 1 

Qi(t)  >  0,  i,t  =  l,...,L. 

Note  that  the  equality  constraints  in  (42)  follow  from 

i  E tr  {rr1/2Qi(*)rr1/2}  =  I E {vf 

^  t= 1  ^  t= 1 

=  zEtr{p^)>  =  p- 

Let  us  denote  Q;(f)  =  UQi(i)Zi(i)U§.(i)  as  the 

eigendecomposition  of  Q^z),  i,t  =  1  where 

Uq; (t)  is  a  unitary  matrix  of  eigenvectors,  and  Z; (z)  = 
dia,g{zii(t),Zi2(t)....,ziN(t)}  \sm  N  x  N  diagonal  matrix 
containing  all  the  eigenvalues.  Then,  the  following  theorem  is 
in  order. 

Theorem  1:  The  optimal  solution  QJ(1), Q*l(L)  of  the 
problem  (41)— (43)  are  diagonal  matrices  and  given  by  Q*(t)  = 
Z*(fUi  =  1, . ,L. 

Proof:  First,  we  note  that  since  S.jj,  i.j  =  1., _ /.  have 

i.i.d.  Gaussian  distributed  entries  with  zero  mean  and  unit  vari¬ 
ance,  the  distribution  of  S ii3-  is  invariant  under  a  unitary  trans¬ 
form  Sv;i?Uq((Z).  Thus,  we  have 


(41) 

(42) 

(43) 


Ia(  Q)  =  7a,(z)  (44) 


where 


s  A  TsT 


Z  =  Z 


FiTl 


1  5  ' 


z;  =  [zf(l),...,zf(P)]T 

A 


Z  i(t)  =  [zii(t)i»i2(t), . . .  rziN(t)Y  i,t  =  l...,L. 


Second,  it  has  been  shown  in  [17]  using  the  theory  of  ma- 
jorization  [18]  that  if  Q ft),  i,t  =  1, ... ,,L  satisfy  the  con¬ 
straint  (42),  then  Z i(t)  also  satisfy  the  constraint 

r{rT1/2zi(i)r?1/2}=p5 

^  t= 1 

Z,(/)  >0:  =  (45) 
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From  (44)  and  (45),  it  can  be  seen  that  Q *(t)  =  Z|(i),  i,t  = 
1, . . .  ,L.  □ 

Using  Theorem  1,  we  can  see  that  since  Q*(t)  = 

Uq.(i)ZJ(i)(UQ.(i))H  =  Z |(i)  for  any  diagonal  matrix 
Z f{t),  obviously,  we  haveUq.(i)  =  Ijy.- 
Substituting  Q f(t)  =  Z*(t)  back  into  (36),  we  obtain  the 
optimal  P i(t)  as 

p* (*)  =  v;r-1/2z*(i)r,-1/2vf ,  i,t  =  i,...,L.  (46) 

We  can  see  from  (46)  that  the  eigenspace  of  the  optimal  source 
covariance  matrix  of  each  transmitting  node  is  equal  to  that  of 
the  corresponding  channel  covariance  matrix.  In  other  words, 
the  channel  should  be  “whitened”  by  the  source  covariance  ma¬ 
trices.  This  property  greatly  reduces  the  number  of  optimization 
variables  from  N2L2  to  NL2.  Now  the  remaining  task  is  to  find 
the  optimal  z.  Itcan  beformulated  as  thefollowing  constrained 
optimization  problem: 


max  Ta(z) 

Z 

(47) 

,  L 

s.t.  jY^^\T-y2Zi(t)T-1/2\  =  Pi 

(48) 

t= 1 

z i(t)  >0,  L  /-T....,  /.. 

(49) 

The  objective  function  (47)  can  be  further  written  as 

l  l  / 

«*)=  EE  Es  {log2  |T i{t)  +  /3i;iSi;iZi(i)S^|} 

^  t=  1  2=1  \ 

— 77s{l°g2  |Tl(7)|}^  .  (50) 

Comparing  (50)  with  (18),  wefind  that  they  have  the  same  struc¬ 
ture.  Thus  (50)  can  be  written  in  thesameclosed  form  as  (20)  by 
following  the  computation  steps  listed  in  Section  III.  Since  the 
constraints  (48)— (49)  are  simple  linear  constraints,  the  problem 
(47)— (49)  can  be  solved  by  the  projected  gradient  technique  de¬ 
veloped  in  Section  IV,  and  hence  more  details  are  omitted. 

VI.  Numerical  Examples 

We  now  illustrate  the  performance  of  the  space-time  power 
scheduling  scheme  presented  earlier.  Forcomparison  in  thecase 
of  i.i.d.  channel,  wewill  consider  thefollowing  five  schemes. 

•  Scheme  A 1:  Space-timeoptimal  power  schedule,  which  is 
based  on  (21)  and  (22). 

•  Scheme  A  2:  Orthogonal  isotropic  power  schedule,  where 
all  links  conduct  transmissions  in  orthogonal  channels 
(such  as  different  times  and/or  different  frequencies),  and 
the  source  covariance  matrix  used  by  all  links  is  propor¬ 
tional  to  the  identity  matrix. 

•  Scheme  A  3:  Uniform  isotropic  power  schedule,  where  all 
links  conduct  transmissions  at  the  same  time  and  the  same 
frequency,  and  the  source  covariance  matrix  used  by  all 
links  is  proportional  to  the  identity  matrix. 

•  Scheme  A4:  Space-only  optimal  power  schedule,  which  is 
based  on  (31)  and  (32). 

•  SchemeA5:  Low  rank  power  schedule  [6],  where  each  link 
uses  a  low  rank  source  covariance  matrix  where  the  corre¬ 
sponding  ranks  for  L  links  are  denoted  by  the  string  of  in¬ 
tegers  (n,r2, . . .  ,rL ).  The  7th  link  with  the  rank  n  uses 


Fig.  5.  Ergodic  capacities  of  schemes  A 1  and  A 4  versus  the  number  of  initial¬ 
izations.  Here,  Ar  =  2,  Z  =  2,  SNR  =  20  dB,  IN  R  =  30  dB. 


a  power  vector  of  n  nonzero  equal  entries  and  N  -  n 
zero  entries. 

For  the  correlated  channel  case,  we  compare  the  following  four 
schemes. 

•  Scheme  Bl:  Space-timeoptimal  power  schedule,  which  is 
based  on  (46)— (49). 

•  Schemes  B  2-B  4:  T  hese  approaches  are  si  mi  lar  to  Schemes 
A 2,  A3,  and  A 5,  respectively.  Flowever,  for  a  fair  compar¬ 
ison,  in  schemes  B2,  B3,  and  B4,  the  design  of  the  source 
covariance  matrices  takes  into  account  the  channel  covari¬ 
ance  matrices. 

Schemes  Al,  A 2,  Bl,  and  B2  are  space-time  based,  and  all 
other  schemes  are  space-only  based.  Schemes  Al  and  Bl  are 
space-timeoptimal,  whileSchemesA2  and  B2  are  not. 

For  Schemes  A 1,  A  4,  and  Bl,  thefollowing  parameters  were 
applied:  s  =  1,  a  =  O.l,  0  =  0.5,  and  e  =  10-4. 

For  all  examples,  we  consider  a  symmetric  and  circular  net¬ 
work  where  all  nodes  are  uniformly  located  on  a  circle.  Fig.  1 
shows  an  example  of  such  a  network  with  L  =  3  active  links. 
For  each  case  to  be  considered,  we  set  the  same  transmit  power 
for  all  transmitting  nodes  P;  =  P,  i  =  1, . . .  ,L,  the  same  gain 
for  all  desired  links  0iyi  =  p,  i  =  1, . . .  ,L,  and  the  same  gain 
for  all  interfering  links  .7,-., •  =  0,  i,j  =  1  ^  /.Thus, 

we  have  the  same  averaged  SNR  for  all  desired  links  Ai  =  Pp, 
i  =  and  the  same  nominal  IN  R  for  all  interfering  links 

Bij  =  P/3,  i,  j  =  1,. . .  ,L,  j  ^  i.  For  all  simulations  of  er¬ 
godic  capacities  versus  I  NR,  we  set  P  =  20  dB  and  p  =  l.The 
INR  is  changed  by  varying  (3.  Similarly,  for  all  simulations  of 
ergodic  capacities  versus  SNR,  we  take  P=  10  dB  and  0=1. 
Then  we  change  p  to  vary  SN  R.  Si  nee  the  nominal  IN  R  and  the 
averaged  SNR  are  quasi  static,  we  assume  that  their  values  can 
beobtained  with  relatively  high  precision.  The  ergodic  capacity 
shown  in  all  figures  is  a  per-link  ergodic  capacity. 

For  Schemes  Al,  A4,  and  Bl,  the  power  allocation  vectors 
were  initialized  randomly.  For  each  simulation  point,  a  certain 
number  of  initializations  were  tried  and  the  best  result  was 
chosen.  Fig.  5  shows  the  ergodic  capacities  of  Schemes  A 1  and 
A4  for  the  i.i.d.  channel  with  N  =  2,  L  =  2,  SNR  =  20  dB, 
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TABLE  II 

Average  Number  of  Iterations;  N  =  2,  L  =  2,  SNR  =  20  dB 


INR  (dB) 

-5 

O 

5 

10 

15 

20 

25 

30 

35 

Scheme  Al 

7 

11 

20 

187 

10 

10 

11 

14 

15 

Scheme  A4 

6 

8 

14 

87 

14 

15 

17 

18 

20 

Fig.  6.  Comparison  of  ergodic  capacities  of  schemes  A1-A5;  i.i.d.  channel. 
Scheme  A 5  uses  (ri,r2)  =  (1,1).  Here,  N  =  2,  L  =  2. 


I N  R  =30  dB,  and  different  number  of  initial  points.  From  this 
figure,  one  can  see  that  the  difference  caused  by  the  number  of 
initializations  is  not  very  large.  In  fact,  after  30  initializations, 
the  difference  is  almost  zero.  Thus,  in  the  sequel,  we  use  30 
initializations  for  Schemes  A 1,  A 4,  and  Bl. 

The  number  of  iterations  required  by  Schemes  Al,  A 4  and 
B 1  for  convergence  depends  on  many  factors  such  as  N,  L, 
SNR,  INR,  and  the  choice  of  the  initial  point.  As  an  example, 
in  Table  II  we  list  the  average  number  of  iterations  for  conver¬ 
gence  of  Schemes  A 1  and  A  4  over  30  random  initializations  for 
different  INR  levels  at  N  =  2,  L  =  2,  and  SNR  =  20  dB.  It  is 
observed  that  the  number  of  iterations  required  at  IN  R  =  10  dB 
is  relatively  large.  It  must  be  due  to  the  insensitivity  of  the  ob¬ 
jective  function  around  the  optimal  covariance  matrices  at  this 
particular  INR  value.  This  INR  value  also  coincides  with  the 
sharp  transition  point  of  the  sum  capacity  of  Scheme  A 1  (see 
Fig.  6). 

I  n  general ,  the  number  of  iterations  for  convergence  increases 
with  the  number  of  optimization  variables  (i  .e.,  NL2).  Flowever, 
an  exact  relationship  of  the  number  of  iterations  with  respect  to 
the  number  of  links  is  difficult  to  obtain. 

The  performance  comparison  of  Schemes  A 1-A  5  for  the  i.i.d. 
channel  environment  are  shown  in  Figs.  6-10,  while  the  perfor¬ 
mance  of  Schemes  B1-B4  for  the  correlated  channel  scenario 
are  compared  in  Figs.  11-15.  Fig.  6  compares  the  ergodic  ca¬ 
pacities  of  schemes  A  1-A 5  with  L  =  2  and  N  =  2.  From 
Fig.  6,  wecan  see  thatSchemeA2  performs  the  sameasScheme 
A 1  at  high  IN  R,  and  Scheme  A3  has  the  same  performance  as 
SchemeAl  atlow  INR.  M  ore  interestingly,  the  transition  of  the 
optimality  from  Scheme  A  2  to  Scheme  A3  along  the  INR  axis 
is  very  sharp  (within  about  1.5  dB  of  INR).  This  property  of 


Fig.  7.  Comparison  of  ergodic  capacities  of  schemes  A1-A5;  i.i.d.  channel. 
Scheme  A  5  uses  (ri,  r2)  =  (1,1).  Here,  N  =  2,  L  =  2. 


Fig.  8,  Comparison  of  ergodic  capacities  of  schemes  A1-A5;  i.i.d.  channel. 
Scheme  A  5  uses  (ri,  r2)  =  (1,1).  Here,  N  =  3,  L  =  2. 


Schemes  A2  and  A3  is  also  observed  in  Figs.  8  and  9.  The  er¬ 
godic  capacities  versus  SN  R  for  the  same  simulation  setups  at 
INR  =  10  dB  is  shown  in  Fig.  7.  Similarly,  we  observe  the  opti¬ 
mality  of  SchemeAl  in  the  whole  SNR  region. 

Fig.  8  compares  the  ergodic  capacities  of  Schemes  A  1-A 5 
withL  =  2  and  TV  =  3.  FI  ere,  Scheme  A  5  with  (n,r2)  =  (1,1) 
remains  strongly  suboptimal  even  compared  with  Scheme  A4 
over  the  whole  range  of  INR.  This  is  because  under  N  =  3, 
there  are  effectively  three  i  ndependent  streams.  B  ut  Scheme  A  5 
with  (n,  r2)  =  (1, 1)  uses  only  two. 

Fig.  9  compares  the  ergodic  capacities  of  Schemes  A  1-A 5 
with  L  =  3  and  TV  =  2.  FI  ere,  the  capacity  of  Scheme  A  5 
with  (ri,r2,r3)  =  (1,1,1)  decreases  as  I N  R  increases.  This 
is  because  it  over  uses  the  total  available  (two)  streams  at  high 
INR. 

Fig.  10  compares  Schemes  A3,  A4,  and  A5,  which  are  all 
space-only  schemes.  FI  ere,  N  =  3  and  L  =  2.  Scheme  A  5 
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INR(dB) 


Fig.  9.  Comparison  of  ergodic  capacities  of  schemes  A1-A5;  i.i.d.  channel. 
Scheme  A  5  uses  (ri,  r2,  r3)  =  (1,1,1).  Here,  N  =  2,  L  =  3. 


INR  (dB) 


Fig.  10.  Comparison  of  ergodic  capacities  of  Schemes  A3,  A4,  and  A5;  i.i.d. 
channel.  Scheme  A  5  uses  (n,  r2)  =  (1, 1),  (ri,  r2)  =  (1, 2),  or  (n,  r2 )  = 
(2,2).  Here,  N  =  3,  L  =  2. 


INR  (dB) 


Fig.  11.  Comparison  of  ergodic  capacity  of  schemes  B1-B4;  correlated 
channel  with  w  =  (0.8, 0.9).  Scheme  B4  uses  (ri,  r2)  =  (1,1).  Here, 
N  =  2,  L  =  2. 


SNR  (dB) 


with  different  rank  conditions  is  illustrated.  Given  N  =  3  and 
L  =  2,  Scheme  A  5  with  (ri,r2)  =  (1,2)  performs  the  same 
as  Scheme  A4  at  high  INR.  When  Ya  ri  <  N,  Scheme  A5 
remains  strongly  suboptimal  compared  with  Scheme  A 4  at  all 
INR. 

Recall  the  optimality  property  of  Schemes  A 2  and  A3  that 
Scheme  A  2  performs  the  same  as  Scheme  A 1  when  I N  R  is  larger 
than  a  threshold  and  Scheme  A  3  has  the  same  performance  as 
Scheme  A1  when  INR  is  less  than  the  threshold.  We  can  de¬ 
termine  the  threshold  INR  value  B*  by  solving  the  following 
nonlinear  equation: 


■7'(f4'iv'£)  = 


J2(j±,N,N) 


which  is  explained  next. 

I  n  (51),  A  =  Pp,  J-iiypp-,  q )  denotes  the  ergodic  sum  capacity 
ofapxr/MIMOlink  with  SNR  v.  Hence,  J2(AL/N,N,N)/L 


Fig.  12.  Comparison  of  ergodic  capacity  of  schemes  B1-B4;  correlated 
channel  with  lj  =  (0.8, 0.9).  Scheme  B4  uses  (r1,r2)  =  (1,1).  Here, 
N  =  2,  L  =  2. 


is  the  ergodic  sum  capacity  using  Scheme  A  2.  U  sing  the  result 
in  [19],  J2(7y,|5,(/)  can  be  written  as 

■h(v,  p,  q) 


=  log2(e)e1/!' 


(—1)^(2 jy.(q-p  +  l)l 
2 2t~ljW.(q  -  p  +  j)\ 


f  2i  —  2j 
\  i~j 


2  j  +  2q-2p 

2  j-l 


q—p+l 

lE^+i 

k= 0 


Here,  Sn(x)  =  e~xrT~ndr,  n  =  0,1,2,...,  is  the  expo¬ 
nential  integral  function  of  order  n  [12,  eq.  5.1.4]  and  can  be 
calculated  using  its  recurrence  relation 


_  e  x  -  xSn-i(x) 
n—  1 
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INR  (dB) 


Fig.  13.  Comparison  of  ergodic  capacity  of  schemes  B1-B4;  correlated 
channel  with  ui  =  (0. 4, 0.5).  Scheme  B4  uses  (ri,r2)  =  (1,1).  Here, 
N  =  2,  L  =  2. 


INR  (dB) 

Fig.  14.  Comparison  of  ergodic  capacity  of  schemes  B1-B4;  correlated 
channel  withu;  =  (0.7, 0.8, 0.9).  Scheme  B4  uses  ( r x ,  r2)  =  (1, 1, 1).  Here, 
jy  =  2,  L  =  3. 


Also,  p  =  mm{NT,NR}  and  q  =  max{NT,NR},  where  NT 
and  Nr  are  the  number  of  transmit  and  receive  antennas,  re¬ 
spectively.  Thus,  p  <  q. 

I  n  (51),  Ji(B*/N,  A/N,  N,  L)  is  the  ergodic  sum  capacity  of 
Scheme  A3.  A  complete  form  of  B (k,  v,  N,  L)  is  as  follows: 

N~  1  1 

=  log2(e)|A|  1  N  A1:Ar(L_i)+i 
2—0  *' 

x£;Ajv(l— i)+i+2:jv_L  —  J-iiy,  N,  NL  -  N ) 

where  A  is  a  matrix  whose  (p,  g)th  element  is  given  by 


AP;9  — 


(^-!)[p  1]  p  =  l,...,N 
(k®-1)  b_Ar_1]  p  =  N  +  1, . . . ,  NL. 


is  a  column  vector  whose pth  entry  is  given  by  (52),  shown  at 
the  bottom  of  the  page,  and  A i:j  stands  for  the  7th  to  jth  col  umns 
of  matrix  A.  H  ere  (•) M  denotes  the  derivative  of  order  i. 

Since  g(B)  =  J1{B/N,A/N1N1L)-J2(AL/N,N,N)/L 
is  a  monotonically  decreasing  function  of  B,  thus  (51)  can 
be  easily  solved,  for  instance,  by  the  bisection  method  [14, 
pp.  145],  As  an  example,  we  applied  (51)  to  the  case  of  Fig.  9, 
and  obtained  B*  =  9.78  dB,  which  agrees  well  with  the 
simulation  result  in  Fig.  9. 

In  practice,  the  threshold  INR  B*  can  be  tabulated  for 
different  network  parameters  such  as  the  number  of  links,  the 
number  of  antennas  of  each  node,  and  the  SNR.  Once  this 
table  is  available,  it  can  be  looked  up  in  real  time  to  determine 


whether  each  node  should  be  scheduled  under  Scheme  A 2  or 
Scheme  A3. 

In  the  following,  we  simulate  the  environment  with  corre¬ 
lated  MIMO  channels  using  the  exponential  correlation  model 
(40).  Fig.  11  compares  the  ergodic  capacities  of  Schemes 
B1-B4  with  L  =  2  and  N  =  2.  FI  ere,  the  channel  is  set  to 
be  highly  correlated  with  u  =  (0.8, 0.9).  For  Scheme  B4, 
we  have  (ri,r2)  =  (1,1).  it  can  be  seen  from  Fig.  11  that 
Scheme  B4  outperforms  Scheme  B2.  Note  that  Scheme  B2 
is  equivalent  to  TDM  A .  We  now  see  that  for  correlated  chan¬ 
nels,  TDMA  is  no  longer  optimal  at  high  INR.  That  is,  the 
channel  correlation  can  change  the  optimality  of  a  particular 
power  schedule  scheme  (exact  for  the  space-time  optimal 
schemes  A 1  and  B 1).  It  can  also  be  observed  that  Scheme  B1 
consistently  provides  the  optimal  ergodic  capacity  throughout 
the  whole  INR  region.  The  optimality  of  Scheme  B 1  is  also 
observed  in  Figs.  13-15. 

The  ergodic  capacities  versus  SN  R  for  the  same  simulation 
setups  at  INR  =  10  dB  is  shown  in  Fig.  12.  Similarly,  we  ob¬ 
serve  the  optimality  of  Scheme  B 1  in  the  whole  SN  R  region. 

Fig.  13  illustrates  the  performance  of  Schemes  B 1-B 4  with 
L  =  2,  N  =  2,  and  low  channel  correlation  u  =  (0.4, 0.5). 
Comparing  Fig.  13  with  Fig.  11,  it  can  be  seen  that  the 
ergodic  capacity  of  Scheme  B2  increases  when  the  channel 
is  less  correlated,  while  the  performance  of  Scheme  B4 
with  (ri,r2)  =  (1,1)  degrades.  The  reason  is  that  when  the 
channel  is  less  correlated,  F,  tends  to  be  an  identity  matrix. 
Thus,  Scheme  B2  is  close  to  be  optimal  in  high  interference 


\  (yN(L  1  jg00  ln(l  +  x)xle  x!v )b  ^  p—  1 , . . . ,  A 

y  (,..v(/.-i)-i  in(!  +  X)xie-X/K)  [p~N~1]  p  =  n  +  1, , . . ,  NL 


(Zi)P  = 


(52) 
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region.  And  Scheme  B4  becomes  suboptimal  when  the 
channel  correlation  is  low.  Interestingly,  for  lj  =  (0.4, 0.5), 
Schemes  B2  and  B4  happen  to  yield  a  similar  ergodic 
capacity  in  the  high  interference  region. 

Fig.  14  compares  the  ergodic  capacities  of  Schemes  B1-B4 
with  L  =  3,  ZV  =  2,  and  a;  =  (0.7, 0.8, 0.9).  We  find  that  in  the 
presence  of  strong  interference,  Scheme  B2  performs  the  same 
as  Scheme  Bl.  The  capacity  of  Scheme  B4  with  (n ,  r>. r3)  = 
(1,1,1)  decreases  as  I  NR  increases,  because  it  over  uses  the 
total  available  streams  at  high  IN  R. 


Finally,  Fig.  15  shows  the  ergodic  capacities  of  Schemes 
B1-B4  with  L  =  2,  ZV  =  3,  and  u  =  (0.8, 0.9).  Interestingly, 
we  can  observe  that  Scheme  B4  consistently  outperforms 
Scheme  B2.  Here,  the  optimality  result  is  different  from 
that  of  the  i.i.d.  channel  environment,  where  the  orthogonal 
isotropic  power  schedule  outperforms  the  low  rank  power 
schedule  (see  Fig.  8).  In  contrast  to  Fig.  11,  Scheme  B4  is 
suboptimal  even  in  the  high  interference  region.  This  fact 
demonstrates  the  importance  of  Scheme  Bl.  The  optimal 
ergodic  capacity  can  not  be  achieved  by  a  “prefixed"  power 


dla(d) 

ddim{t) 


=  1°^  E  E  «iknQ(n,  a t,i,k)  +  <hiknQ'(n,  AUfe))  +  1°^  E 
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Fig.  15.  Comparison  of  ergodic  capacity  of  schemes  B1-B4;  correlated 
channel  with  u  =  (0.8, 0.9).  Scheme  B4  uses  (r!,r2)  =  (1,1).  Here, 
N  -  3,  /.  -  2. 


scheduling  schemes.  Our  schemes  A1  and  B1  are  adaptive  to 
the  network  conditions. 


VII.  Conclusion 

We  have  proposed  a  space-time  optimal  power  scheduling 
approach  for  multiple  distributed  MIMO  links  assuming  no 
instantaneous  CSI  at  the  transmitting  nodes.  For  the  i.i.d.  and 
correlated  channel  environments,  thisapproach  leads  to  Scheme 
A 1  and  Scheme  B 1,  respectively.  I  n  particular,  for  a  symmetric 
network  of  i.i.d.  channels,  using  Scheme  A1  as  the  optimal 
benchmark,  we  have  observed  that  Scheme  A  2,  an  orthogonal 
isotropic  power  schedule  (such  as  TDM  A),  is  optimal  when 
the  I N  R  is  larger  than  a  threshold,  and  Scheme  A  3,  a  uniform 
isotropic  power  schedule,  is  optimal  when  the  IN  R  is  less  than 
the  same  threshold.  The  threshold  IN  R  value  can  be  computed 
based  on  the  network  topology  and  the  transmit  power  level. 
For  the  correlated  channel  environment,  we  have  shown  that  the 
channel  correlation  has  a  great  impact  on  the  ergodic  capacity 
and  the  optimality  of  different  power  scheduling  approaches. 
Future  research  should  address  issues  arising  from  asymmetric 
networks. 


Appendix  A 

Summary  of  Notations  for  System  Model 

H  ij  MIMO  channel  between  the  2th  transmitting  node 

and  the  jth  receiving  node. 
i  =  j,  channel  gain  of  a  link-of-interest;  i  ^  j, 
channel  gain  of  an  interfering  link. 

Pi(i)  Transmit  covariance  matrix  of  the  «th  transmitting 
node  in  the  vth  time  slot. 

R;(i)  Interference-plus-noise  covariance  matrix  at  the* 
receiving  node  in  the  tt h  time  slot. 

Pi  Averaged  transmit  power  from  the  ith  transmitting 
node. 


Xi(i)  Signal  vector  transmitted  from  the  «th  transmitting 
node  within  time  slot  t. 

y i(t.)  Signal  vector  received  by  the-ith  destination  node 

within  time  slot  t. 

Ia  Averaged  ergodic  sum  capacity. 

I  Ergodic  sum  capacity  at  one  time  slot. 

fli  Channel  correlation  matrix  of  the-ith  transmitting 
node. 


Appendix  B 
The  Gradient  of  (20) 

See  the  first  equation  shown  at  the  bottom  of  the  previous 
page,  where  (•)'  denotes  _the  first  order  derivative,  dtikn,  u[ikn, 
Q'(n,Xt,i,k),  and  Q'(n,Xtiitk )  are  given  by  (53)-(56),  respec¬ 
tively,  shown  at  the  bottom  of  the  previous  page.  FI  ere,  q  = 
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